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The system size dependence of the multifractal spectrum f{a) and its singularity strength a 
is investigated numerically. We focus on one-dimensional (ID) and 2D disordered systems with 
long-range random hopping amplitudes in both the strong and the weak disorder regime. At the 
macroscopic limit, it is shown that /(a) is parabolic in the weak disorder regime. In the case of 
strong disorder, on the other hand, /(a) strongly deviates from parabolicity. Within our numerical 
uncertainties it has been found that all corrections to the parabolic form vanish at some finite value 
of the coupling strength. 
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I. INTRODUCTION 

In recent years, it has become clear that the eigenfunc- 
tions of noninteracting disordered systems at the critical 
point of a metal-insulator transition (MIT) have multi- 
fractal scaling properties. Hence, multifractal analysis 
has become a standard tool for treating the strong wave- 
function fluctuations which arise near criticalityji One 
way to characterize wave-function statistics is through 
the corresponding multifractal /(a) spectrum. This spec- 
trum of the local density is universal at the critical point, 
and the exponents such as ao, which describes the scal- 
ing behavior of the typical local electron density, are suit- 
able for describing the phase transitioniS Finding the true 
/(a) curve is essential in order to fully understand MITs. 

MITs depend on the dimensionality and symmetries 
of the system and can occur in both the strong disor- 
der and the weak disorder regime (strong-coupling or 
weak-coupling regime, respectively, in the language of 
field theory). Each regime is characterized by the cor- 
responding coupling strength that depends on the ratio 
between diagonal disorder and the off-diagonal transition 
matrix elements of the Hamiltonian.^ We wish to point 
out that most of the already published analytical results 
concerning the f{a) spectrum, which we summarize in 
the next paragraphs, have been obtained in the weak 
disorder limit (weak-coupling approach). 

Using the renormalization group in d = 2 + e dimen- 
sions {e = ao — d <^ 1), it was found^ that the singularity 
spectrum is parabolic, 



f{a)=d 



{a - gp)^ 
4(ao - d) 



(1) 



up to O(e^). This equation was justified for the case of 
unbroken time-reversal symmetry (orthogonal universal- 
ity class) in the weak-coupling limit. 

Recently, in Ref. 5 it was found that a parabolic spec- 
trum (PS) is fulfilled for the critical wave function of a 
two-dimensional (2D) Dirac fermion in a random vector 
potential (chiral universality class). These results were 
obtained by mapping the problem to a Gaussian field the- 
ory in an ultrametric space (a Cayley tree). For a deep 



understanding of the origin of the multifractal behavior 
in this model, see also Ref. 6. 

The field-theoretical approaches in the context of the 
integer quantum Hall transition (unitary universality 
class) has been addressed more recently. Based on the 
application of conformal field theory to the critical point 
of a 2D Euclidean field theory, Bhaseen et al^ conjec- 
tured that the singularity spectrum is exactly parabolic 
at the plateau transitions. Using the same general prin- 
ciples of conformal field theory, '^•^ a PS was conjectured 
for a closely related quantity, the two-point conductance. 

Given the present status of this problem, it would be 
highly desirable to check (i) whether the PS is fulfilled 
at weak disorder strength and (ii) whether or not the PS 
holds in the strong disorder regime, where it is known 
that all realistic MITs, such as the conventional Ander- 
son transition in 3D, take place. Unfortunately, because 
of the absence of a small parameter, no quantitative pre- 
dictions for the critical properties and, in particular, for 
the / (a) spectrum have been made in the strong-coupling 
regime. In the absence of such predictions, the problem 
must be addressed using numerical calculations. That is 
why we are interested in exploring the role of the disor- 
der strength on the /(a) spectrum. Specifically, we want 
to analyze how the singularity spectrum evolves in the 
face of increasing coupling strength, i.e., when the tran- 
sition shifts from the weak disorder regime to the strong 
disorder regime. 

This paper is organized as follows. In Sec. I, we in- 
troduce the model and the methods used for the calcu- 
lations. The results for the multifractal spectrum in ID 
and 2D are shown in Sees. II A and II B, respectively. 
Finally, the conclusions are presented in Sec. HI. 



II. MODEL AND METHODS 

For computation of the multifractal spectrum, we use 
the standard box-counting procedure^ first dividing the 
system of i"* sites into Ni ~ {L/lY boxes of linear size I 
and determining the box probability of the wave function 
in the i box by = X]„ IV'fenPj where the summation 
is restricted to sites within that box and ipun denotes the 



2 



amplitude of an eigenstate with energy at site n. The 
normaUzed qth moments of this probabiUty fj,i(q,l) — 

A*i (0/ Sjlli A*j(0 constitute a measure. From this, the 
Lipschitz-Holdcr exponent or singularity strength can be 
obtained 

«^(L).lim^£ii^iM^^^Mili), (2) 
5^0 Ino 

as well as the corresponding fractal dimension 

f{ag{L)) = hni t — , (3 

A^o Ino 

which yields the characteristic singularity spectrum /(a) 
in a parametric representation. In Eqs. l|2Jl and 
S = l/L denotes the ratio of the box sizes and the sys- 
tem size. If the qth moments of the measure counted 
in all boxes are proportional to a power Tq of the box 

size, (X^il^i A^i (0) °^ l"^"? multifractal behavior may be 
derived. The /(a) spectrum and the mass exponent 
Tq are related by a symmetric Legendre transformation, 
/(a) ^ aq ~ Tq, with a = dTq/dq and q ~ df{a)/da. 
From Eq. (Q, one can easily obtain 

aQ — d = d — ai . (4) 

In what follows, we will first use this relation to check 
whether or not the PS is valid in each regime stud- 
ied. First, aq{L) and f{aq{L)) were calculated for 
different system sizes and then extrapolated to the 
macroscopic limit aq = lim^^oo aij(i) and f{aq) — 
linii^oo f{aq{L)). We wish to clarify that the calculation 
of aq{L) and f{aq{L)) is suitable only if the conditions? 

a < / < L < ^ (5) 

are satisfied, where ^ is the localization or correlation 
length and a is the lattice spacing (or any microscopic 
length scale of the system, such as the cyclotron radius) . 

Usually, the disorder-induced MIT is investigated for 
Hamiltonians with short-range, off-diagonal matrix ele- 
ments (e.g., the canonical Anderson model). Another 
class of Hamiltonian exhibiting an MIT in arbitrary di- 
mension d is formed by those that include long-range 
hopping terms. The effect of long-range hopping on lo- 
calization was originally considered by Andersor>ii for 
randomly distributed impurities in d dimensions with 
the V{r — r') ~ |r — r'\~^ hopping interaction. It 
is knowniiii^ that all states are extended for (3 < d, 
whereas for (3 > d, the states are localized. Thus, one 
can study the MIT by varying the exponent /? at fixed 
disorder strength. At the transition line (3 — d, & real- 
space renormalization group can be constructed for the 
distribution of couplings (iSiH These models are most con- 
venient for studying critical properties because the exact 
critical point is known (/3 = d) and they allows us to treat 
the ID and 2D cases, thus reaching larger system sizes 



and reducing the numerical effort. Here, we will concen- 
trate on the ID and 2D versions of these models with 
orthogonal symmetry. In our calculations, we consider a 
small energy window, containing about 8% of the states 
around the center of the spectral band. The number of 
random realizations is such that the number of critical 
states included for each L is roughly 4 x 10^, while in or- 
der to reduce the edge effects, periodic boundary condi- 
tions in all directions are included. Using methods based 
on level statistics, we checked that the normalized near- 
est level variancesi^ are indeed scale invariant at each 
critical point studied. 

A. ID system 

First, we concentrate on the intensively studied power- 
law random banded matrix model (PRBM)iS*iSiiLi2iia 

The corresponding Hamiltonian that describes a disor- 
dered ID sample with random long-range hopping is rep- 
resented by real symmetric matrices, whose entries are 
randomly drawn from a normal distribution with zero 
mean, {TLij) = 0, and a variance which depends on the 
distance between lattice sites, 

The model describes a whole family of critical theories 
parametrized by < 6 < oo, which determines the crit- 
ical dimensionless conductance in the same way as the 
dimensionality labels the different Anderson transitions. 
In the two limiting cases 6^1 and 6 <C 1, which 
correspond to the weak and the strong disorder limits, 
respectively, some critical properties have been derived 
analyticallyii^*iSiS242ii2^ The system size ranges between 
L = 40 and L = 1800, and 0.03 < 6 < 30. We restrict 
ourselves to values oi q > 0, which correspond to a < ctQ. 

Using the exact eigenstates of Eq. © from numeri- 
cal diagonalizations, we obtained the size dependence of 
the scaling exponents OLq[L) and f{aq{L)) at the criti- 
cal point of the finite system. In order to extrapolate to 
macroscopic systems, we propose a finite-size correction 
to aq{L) = Uq + aqL^^i , where the irrelevant exponent 
yq > 0. For all values of q, we found that the exponent 
Tjq hardly differs from unity. The same behavior is found 
for f[aq{L)). Thus, we can write 

aq{L) = aq + aq/L , f{aq{L)) = f{aq) + bq/L, (7) 

with aq, aq, f{aq), and bq being adjustable parameters. 
As we will see below, our numerical data strongly support 
this behavior. 

As a first step, for each value of g, we evaluate the 
numerators on the right-hand sides of Eqs. ||2J) and ||2J), 
respectively, for decreasing box sizes, and we calculate 
aq{L) and f{aq{L)) from the slopes of the graphs of 
the numerators versus In (5. In order to satisfy the in- 
equalities of Eq. ||SJ), we take S to be in the interval 
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FIG. 1: iJ.i{q,l)ln as a function of InS for q = 

(open symbols, right axis) and g = 1 (solid symbols, left axis) 
and different system sizes: L = 40 (circles), 60 (triangles), 
120 (diamonds), and 1800 (squares). The straight lines whose 
slopes correspond to the values of ao{L) and ai(L) are linear 
fits to Eq. iP. 
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FIG. 2: The L dependence of the scaling exponents ao (cir- 
cles) and ai (squares) for the ID disordered system in the 
strong-coupling regime {b — 0.3). The straight lines are lin- 
ear fits to Eq. 0. The inset shows the same dependence for 
the case of weak coupling {b = 10). 



(0.1,0.5). Figure 1 provides an example of the linear 
fit to fJ-iiq, I) liiMiCli versus InS for two values of q, 
q — (open symbols, right axis) and q = I (solid sym- 
bols, left axis), and different system sizes: L = 40 (cir- 
cles), 60 (triangles), 120 (diamonds), and 1800 (squares). 
Clearly, there is no ambiguity in the determination of the 
slopes that correspond to the values of ao{L) and ai{L). 
These slopes, which are summarized in Fig. 2, have been 
obtained for b = 0.3. 

In Fig. 2, we represent the finite-size corrections for 
the scaling exponents ao (circles) and ai (squares) for the 
ID disordered system described by Eq. © in the strong- 
coupling regime (6 — 0.3). The straight lines are linear 
fits to Eq. lO and intercept the vertical axis at 0.590 and 
0.455, respectively. These results clearly demonstrate 
that Eq. Q is not fulfilled at all, indicating that the 
PS is no longer valid in this case. This is not surprising 
since the parameter found (e = 0.590), is not small, and 
the e expansion can only be justified parametrically 
for e ^ 1. We have found that if b is reduced, the in- 
tercept points shift to higher values and their separation 
increases. 

The inset of Fig. 2 shows the same L dependence for 
the case of weak coupling {b = 10). Unlike in the strong- 
coupling regime, both lines intercept the y axis practi- 
cally at the same values, of 0.0194 and 0.0193, respec- 
tively. This suggests that f{a) is parabolic, in agreement 
with Eq. (Q). In this case, e — 0.0194 <^ 1, as expected. 
Very recently, results similar to those in the inset were 
obtained for the Chalker-Coddington model, although in 
this case a different method^'^ was used to evaluate aq{L). 



In this calculation, the parameter found was e = 0.26. 

The extrapolated results of the whole /(a) curve in 
the strong-coupling regime 6 = 0.1 are depicted in Fig. 
3 (solid circles). The implicit parameter q ranges from 
to 1, in steps of 0.1. Note that Eq. ^ is assumed to 
be valid least up to \q\ < 1.^ The PS with ao = 2.228 
is shown for comparison (solid line). The difference be- 
tween both results is noteworthy. Based on a resonance 
approximation, the following formula was derived in Ref. 
16 for f{a) in the regime 6^1 for values oi q > i, 

f{a) = 2bF{a/2b) , (8) 

where F{A) is the Legendre transform of the function 
r{2q - l)/229-3r(g)r(g - 1). The result of Eq. © for 
= 0.1 is represented by a dot-dashed line. At q — ^, 
the mass exponent Tq diverges, while for smaller values 
of (J, this approximation completely breaks down. Note 
that /(a) properly interpolates between the parabola and 
Eq. |(SJ). For the smallest values of q near the maximum, 
/(a) is very close to the PS, and as q increases, f{a) 
separates from this curve with a regular tendency until 
it reaches the resonance prediction for the larger values 
of q as shown. We have found that as b decreases, the 
deviations with respect to PS are larger, and as b in- 
creases, these deviations diminish. For values of b well 
inside the weak-coupling domain (see inset), the devi- 
ations are practically unobservable. The inset of Fig. 
3 shows the same f{a) spectrum for the weak coupling 
regime (6 = 10). The perfect agreement between the ob- 
tained results (solid circles) and the PS, Eq. with 
ao = 1.0194 is evident. Similar results were obtained in 
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FIG. 3: The f{a) spectrum (solid circles) for the ID disor- 
dered system in the strong-coupling regime (6 = 0.1). The im- 
plicit parameter q ranges from to 1. The PS with ao = 2.228 
is shown for comparison (solid line). The dot-dashed line cor- 
responds to Eq. (|HJ for q > 1/2. The inset shows the same 
spectrum in the case of weak-coupling {b = 10) together with 
the PS with ao = 1.0194 (solid line). 



Ref. 16 for this regime. We can thus conclude, as we 
previously stated, that the PS is correct in this case. 



B. 2D system 

We have also calculated the /(a) spectrum in the ex- 
perimentally more important case d = 2. Unlike the 
ID PRBM model, until now, it has not been possible to 
solve the 2D disordered models with long-range transfer 
terms analytically. We consider nonintcracting electrons 
on a 2D square lattice with random on-site potentials Si 
and random transfer terms Vij — Veij/\ri — rj p, whose 
Hamiltonian reads 



ri){ri 



(9) 



where the vectors label the sites of the lattice, and {si] 
and {sy} are two sets of uncorrelated random numbers 
uniformly distributed within the interval (— VF/2,W/2), 
and {~S/2, S/2), respectively. V = 1 defines the energy 
scale and S/V is taken to be equal to 1 in all regimes. 
Each regime is characterized by the coupling strength 
A = W/S and we take 0.4 < W/V < 30. The system 
sizes used are L — 24, 36, 56, 72, and 96. The critical 
properties of this model in the strong disorder regime 
have recently been investigated numerically.^^ A closely 
related model in c? = 1 — 3 was also considered in Ref. 25 
for strong disorder. 



FIG. 4: As for Fig. 2, for the 2D disordered system in 
the strong-coupling regime (A = 6). The results of the inset 
correspond to the same system in the weak-coupling regime 
(A = 0.4). 



In Fig. 4, we show the same size corrections as in Fig. 
2 for the 2D disordered system in the strong-coupling 
regime (A = 6). The fitted straight lines intercept the 
vertical axis at 0.827 and 0.743. The corresponding re- 
sults for the case of weak coupling (A = 0.4) are reported 
in the inset. As in the ID system, we obtained a be- 
havior for these corrections in both regimes. All the com- 
ments made with respect to the ID system are equally 
valid in this case. Hence, we arrive at the same con- 
clusions for the 2D case, i.e., the PS is valid in the weak 
disorder regime whereas the strong disorder case strongly 
deviates from this behavior. 

Figure 5 shows the same extrapolated values of the 
/(a) spectrum as in Fig. 3 for the 2D disordered system 
in the strong disorder regime (A — 30). The PS with ao = 
4.3673 is shown for comparison (solid line). As in the ID 
case, the difference between both results is notable. As 
there are no analytical predictions for this regime, we are 
not able to justify the sign of the deviations in this case. 
The inset shows the same f{a) spectrum (solid circles) 
for the weak-coupling regime (A = 0.4) together with the 
PS (sohd hue) with ao = 2.0692. 

We now define in a simple way the extent to which the 
spectra deviate from parabolicity. We consider the dis- 
tance between the two intercept points in Figs. 2 and 4, 
^ad = ao + ai — 2(i, as a measure of departure from the 
parabolic shape. This clearly distinguishes a parabolic 
spectrum from a nonparabolic one. In the following, we 
will use the notations ti — 1/b and t2 ^ \ for the cou- 
pling constant of the ID and 2D models, respectively. 
In Fig. 6, we represent the dependence with the cou- 
pling constant of the parameter Aad for ID (squares) 
and 2D (circles) disordered systems. This parameter will 
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FIG. 5: Same as Fig. 3, for the 2D disordered system in the 
strong-coupling regime (A = 30). The PS with ao ~ 4.3673 is 
shown for comparison (solid line). The inset shows the same 
spectrum in the case of weak coupling (A = 0.4) together with 
the PS with ao = 2.0692 (solid line). 
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FIG. 6: The small coupling constant dependence of the pa- 
rameter Add for the ID (squares) and 2D (circles) disordered 
systems. The inset shows the same dependence in the whole 
range of 6 and A studied; solid lines are fits to Eq. TM . 



tend to a well-defined finite value in the limit — s- cx). 
The small i^-dependent behavior of Aad, corresponding 
to the weak-coupling regime, is shown in the main panel 
of Fig. 6, in which one can clearly appreciate that below 
a given value of (ti « 0.1 for the ID system and i2 ~ 4 
for the 2D case), A^d = within the error bars. This 
result strongly suggests that all corrections to parabolic- 
ity vanish at a finite value of the coupling constant, in 



contrast to what might be expected naively. 

In order to predict the large asymptotic values of 
Aad from the inset of Fig. 6, a curve of the form 
Aad = Ao — GA,B(td) was fitted to the data points in 
this graph, where Aq, A, and B are three positive fitting 
parameters, and Ga,b {td) is a function which goes to zero 
astd ^ oo. Note that this fit was performed only over the 
larger values of td where convergence was evident. Var- 
ious forms for GA,B{td) were chosen and tested on the 
data plots. These included exponential, inverse logarith- 
mic, and power-law decay with td- The testing involved 
performing a three-parameters fit on the Aad versus 
td data plots using the standard Levenberg-Marquardt 
method for nonlinear fits. The form of GA,B{td) even- 
tually chosen was the exponential law. The reason for 
choosing the above form of GA.Bitd) with preference to 
the alternative forms was simply that the inverse loga- 
rithmic and the power laws do not fit our data properly. 
Thus, the small and large id-dependent behavior of Aq.^ 
can be described by 



A„ 





Ao 



Ae 



h <0.1, i2 <4 

id > 1 ■ 



(10) 



The solid lines in the inset of Fig. 6 are fits to Eq. I|10|) 
with parameters Ao = 1.19 ± 0.02, A = 1.250 ± 0.013, 
and B = 0.054 ± 0.002 for the ID system, and Aq = 
0.88 ±0.02, A = 1.110 ±0.013, and B = 0.055 ±0.002 for 
the 2D case. 

Very recently, Potempa and Schweitzer— reported the 
/(a) spectrum for a critical wave function of a finite sam- 
ple in 2D of linear size L — 150 for A = 6. For this value 
of A, we found Aa2 = 0.084, indicating that we are very 
close to the regime where deviations are practically un- 
observable. These authors found deviations (downwards) 
from PS mainly for the largest negative q values. We be- 
lieve that these deviations are due to the use of the finite- 
size value of OfQ. If the extrapolated value of ao = 2.827 
is used instead, the agreement with the PS is fairly good. 

Before concluding, we must admit that experiments 
to measure ao or related multifractal exponents have 
not been carried out. However, the statistical properties 
of wave functions are close related to the conductance 
statistics. As was demonstrated in Ref. 9, multifractal- 
ity can show up in transport experiments. Specifically, 
the gth moments of the two-point conductance exhibit 
multifractal statistics. The connection between the PS 
for the conductance and Eq. is through the relation 
Xt — 2(q!o — d), where Xt is the power-law exponent of 
the typical conductance. In this way, our results could 
be checked in experiments with random systems using 
samples with varying disorder. 



III. CONCLUSIONS 

In this paper, we have calculated the multifractal scal- 
ing exponents of the wave functions for ID and 2D dis- 
ordered systems with long-range transfer terms at criti- 
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cality. The leading finite-size corrections to a and f{a) 
decay algebraically with exponents equal to —1. In the 
limit L — > cx), we have demonstrate that according to 
theoretical predictions, the /(a) spectrum is parabolic 
in the weak disorder regime. Our calculations strongly 
suggest that the parabolic spectrum completely breaks 
down for large values of the disorder strength, and more 
importantly, that all corrections to the parabolic shape 
vanish (within our numerical uncertainties) at a finite 
value of the coupling constant. 

The question arises as to whether these results are ap- 
plicable to other quantum systems, particularly in the 
3D Anderson transition that occurs in the strong disor- 
der domain, and whose similarity with the PRBM model 
at & = 0.3 has been demonstrated for several critical 



magnitudesji^ Another interesting question is whether 
the deviations studied here also have signatures in other 
properties, such as spectral statistics or generalized di- 
mensions. Finally, how the present 2D model can be 
studied, using 2D conformal field theory or locator ex- 
pansion instead of the resonance pair approximation, re- 
mains to be resolved. 



Acknowledgments 

The author thanks the Spanish DGESIC for financial 
support through Project No. BFM2000-1059. 



* Electronic address: lecr@um.esl 
URL; http : //bohr . f cu.um. es/miembros/ecr/l 

^ V.I. Falko and K.B. Efetov, Europhys. Lett. 32, 627 

(1995); Phys. Rev. B 52, 17413 (1995). 
2 M. Janssen, Int. J. Mod. Phys. B 8, 943 (1994); B. Huck- 

estein. Rev. Mod. Phys. 67, 357 (1995). 
^ K.B. Efetov, Adv. Phys. 32, 53 (1983). 

C. Castellani and L. Peliti, J. Phys. A 19, L429 (1986); F. 

Wegner, Nucl. Phys. B 280, 210 (1987). 
^ C.C. Chamon, C. Mudry, and X.-G. Wen, Phys. Rev. Lett. 

77, 4194 (1996); H.E. Castillo, C.C. Chamon, E. Fradkin, 

P.M. Goldbart, and C. Mudry, Phys. Rev. B 56, 10668 

(1997). 

® David Carpentier and Pierre Le Doussal, Phys. Rev. E 63, 
026110 (2001). 

M.J. Bhaseen, I.I. Kogan, O.A. Soloviev, N. Taniguchi, 
and A.M. Tsvelik, Nucl. Phys. B 580, 688 (2000). 

* M.R. Zirnbauer, hep-th/9905054 (unpublished). 

® M. Janssen, M. Metzler and M.R. Zirnbauer, Phys. Rev. 

B 59, 15836 (1999). 
^'^ Ashvin Chhabra and Roderick V. Jensen, Phys. Rev. Lett. 

62, 1327 (1989). 
" PW. Anderson, Phys. Rev. 109, 1492 (1958). 
^2 L.S. Levitov, Europhys. Lett. 9, 83 (1989); Phys. Rev. 

Lett. 64, 547 (1990). 
" L.S. Levitov, Ann. Phys. (Leipzig) 8, 697 (1999). 



" E. Cuevas, Phys. Rev. Lett. 83, 140 (1999); E. Cuevas, E. 
Louis, and J. A. Verges, ibid. 77, 1970 (1996). 
A.D. Mirlin, Phys. Rep. 326, 259 (2000); F. Evers and 
A.D. Mirlin, Phys. Rev. Lett. 84, 3690 (2000). 

" A.D. Mirlin and F. Evers, Phys. Rev. B 62, 7920 (2000). 

E. Cuevas, V. Gasparian and M. Ortuno, Phys. Rev. Lett. 
87, 056601 (2001); E. Cuevas, Phys. Rev. B 66, 233103 
(2002). 

^® E. Cuevas, M. Ortuno, V. Gasparian, and A. Perez- 

Garrido, Phys. Rev. Lett. 88, 016401 (2002). 

I. Varga and D. Braun, Phys. Rev. B 61, R11859 (2000); 

Imre Varga, ibid. 66, 094201 (2002). 
^° V.E. Kravtsov and K.A. Muttalib, Phys. Rev. Lett. 79, 

1913 (1997). 

A.D. Mirlin, Y.V. Fyodorov, P.M. Dittes, J. Quezada, and 

T.H. Seligman, Phys. Rev. E 54, 3221 (1996). 

V.E. Kravtsov and A.M. Tsvelik, Phys. Rev. B 62, 9888 

(2000). 

F. Evers, A. Mildenberger, and A.D. Mirlin, Phys. Rev. B 
64, 241303(R) (2001). 

H. Potempa and L. Schweitzer, Phys. Rev. B 65, 
201105(R) (2002). 
'^^ D. A. Parshin and H. R. Schober, Phys. Rev. B 57, 10232 
(1998). 



